Selectivity of the ^^0(e,e'pp) reaction to discrete final states 

C. Giusti and F. D. Pacati 
Dipartimento di Fisica Nucleare e Teorica dell'Universita, Pavia 
and Istituto Nazionale di Fisica Nucleare, Sezione di Pavia, Italy 

K. AUaart and W. J. W. Geurts 
Department of Physics and Astronomy, Vrije Universiteit, 
De Boelelaan 1081, 1081 HV Amsterdam, The Netherlands 

W. H. DickhofT 

Department of Physics, Washington University, St. Louis, MO 63130, USA 

H. Miither 

Institut fiir Theoretische Physik, Universitdt Tubingen, 
Auf der Morgenstelle 14, D-72076 Tiibingen, Germany 
(February 9, 2008) 

Resolution of discrete final states in the ^^0(e,e'pp)^*C reaction may provide an interesting tool to 
discriminate between contributions from one- and two-body currents in this reaction. This is based 
on the observation that the 0^ ground state and first 2^ state of ^''C are reached predominantly by 
the removal of a ^ So pair from ^^0 in this reaction, whereas other states mostly arise by the removal 
of a pair. This theoretical prediction has been supported recently by an analysis of the pair 
momentum distribution of the experimental data Q. In this paper we present results of reaction 
calculations performed in a direct knock-out framework where final-state interaction and one- and 
two-body currents are included. The two-nucleon overlap integrals are obtained from a calculation 
of the two-proton spectral function of '^^O and include both long-range and short-range correlations. 
The kinematics chosen in the calculations is relevant for recent experiments at NIKHEF and Mainz. 
We find that the knock-out of a proton pair is largely due to the (two-body) A-current. The 
pair knock-out, on the other hand, is dominated by contributions from the one-body current 
and therefore sensitive to two-body short-range correlations. This opens up good perspectives for 
the study of these correlations in the ^®0(e,e'pp) reaction involving the lowest few states in ^'*C. 
In particular the longitudinal structure function /oo, which might be separated with super-parallel 
kinematics, turns out to be quite sensitive to the NN potential that is adopted in the calculations. 

PACS numbers: 21.10Jx, 21.30.Fe, 21.60.-n, 25.30. Fj. 



I. INTRODUCTION 



Exclusive (e,e'pp) reactions on nuclei have recently been added to the rich set of tools exploring the nucleus with the 
electromagnetic interaction . It is expected that this new tool may finally be able to clarify the nature and influence 
of short-range correlations (SRC) in low-energy nuclear phenomena. Several early theoretical papers [||,|| established 
a link between two-nucleon removal cross sections and the two-nucleon density matrix ||^ or the two-nucleon spectral 
function [pj w hich contain information related to SRC. A somewhat different perspective to this issue has been explored 
in Refs. Mfcl. The anticipated availability of this reaction generated renewed theoretical interest p-pl in the reaction 
description and the calculation of the two-nucleon spectral function. In a series of publications [p|-|l3| a reaction 
description of electron- and photon-induced two-nucleon emission processes has been established. Indeed, it appears 
from these studies that the most promising reaction to study short-range phenomena involves the (e,e'pp) channel, 
where the effect of meson-exchange currents and A-isobars is less dominant as compared to the (e,e'pn) and (7,NN) 
processes. 

In its most recent form the description of the (e,e'pp) excitation process includes the contribution of the usual 
one-body terms as well as those two-body currents which involve the intermediate excitation of the A-isobar. The 
de-excitation of the A after absorption of the photon or the excitation of the A before absorption of the photon 
proceeds by exchange of a pion with another nucleon. In the present work an improvement of the dynamic aspects 
of the propagation of the A-isobar is taken into account iQ. A treatment of A-propagation involving the exchange 
of rho-mesons is not included at present. The treatment of the final state interaction of the outgoing protons with 
the remaining nucleus is treated by neglecting their mutual interaction but including the distorting effect of their 
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interaction with the remaining nucleons in terms of an optical potential. The latter distortion of the individual protons 
is constrained by experimental data obtained from elastic scattering of nucleons off nuclei. The approximation to 
neglect the interaction between the two outgoing protons has been justified in the past by arguing that the pair of 
protons will leave the nucleus largely back to back making this type of final state interaction less important. This 
issue should be further studied in the future since there is no a priori dominance of the effects of correlations before 
or after the absorption of the photon as emphasized in Ref. 1^. It is, however, possible that angular momentum and 
parity restrictions associated with the transition to specific discrete final states in the remaining nucleus, may filter 
the importance of this type of final state interaction. The possibility to analyze different final states in the reaction 
has already been explored in Ref. . The separation of some of the low-lying final states has recently been realized 
experimentally at the NIKHEF and Mainz [g5| facilities. As will be shown in this paper, this feature plays a crucial 
role in clarifying the distinction between transitions that are strongly influenced by SRC and those where two-body 
transition currents play a dominant role. A related description of the reaction process has been presented in Refs. 

.... 

The critical information about SRC in the transition to the final A-2 state is incorporated in the two-body spectral 
function at the corresponding energy. At low missing energy, it represents the probability density for the removal of a 
pair of nucleons (protons in the present work) from the ^^O ground state to a specific discrete final state in -'^^C. Since 
this removal amplitude involves nucleons close to the Fermi energy, the accurate description of this process requires 
a careful treatment of the influence of low-energy, or long-range, correlations associated with the soft-surface features 
of the ^^O nucleus. The latter feature has not been included in Ref. [|l2[, but is incorporated in Ref. It is the 
purpose of the present work to combine the reaction description of the two-proton removal process of Ref. p2| with 



the many-body calculation of the two-particle spectral function in ^^O of Ref . |18| in order to calculate cross sections 
for the triple-coincidence experiments performed at NIKHEF and Mainz. 



The calculation of the two-body spectral function in Ref. |18 includes the dressing of individual nucleons through 
their coupling to low- lying core excitations. In addition, the reduced presence of these nucleons at low energy associated 
with strength removal due to the influence of SRC is incorporated |^ . This yields theoretical spectroscopic factors 
for low-lying states in ^^N which represent the closest agreement with experiment to date. Consistency between 
the two aspects of the calculations (long-range vs. short-range) is ensured by employing the same effective interaction 
(G-matrix) in the calculation of the long-range correlations which is responsible for the removal of single-particle 
strength. Although the appearance of high- momentum nucleons in the ground state is implied by SRC, their presence 
is only apparent at high excitation energy in the A-1 system |^,^. The corresponding cross section for the removal 
of high-momentum protons from ^^O in the (e,e'p) reaction has recently been calculated in Ref. ||2^. Although these 
cross sections are large enough to be detectable at these high energies, other competing processes will also be present 
making a clear-cut identiflcation of SRC in the (e,e'p) reaction difficult. 

This elusive consequence of SRC in the (e,e'p) reaction does not pertain to the removal of a pair of nucleons leading 
to a discrete final state in the A-2 system since few other competing processes are present. The strongly reduced 
probability for a pair of protons to be in close proximity will unavoidably lead to the presence of high-momentum com- 
ponents in their relative momentum wave function. The character and strength of these high-momentum components 
depends on certain aspects of short-range phenomena which are described differently by different nucleon-nucleon 
(NN) interactions. Sensitivity to the choice of the NN interaction in describing pairs with high relative momentum in 
the two-body spectral function has been established in Ref. . It is hoped that a realistic treatment of the reaction 
process combined with a detailed many-body treatment of the spectral function in conjunction with new experimental 
data will produce the first clear and unambiguous determination of SRC in nuclei. 

As will be discussed in this work, the presence of several different discrete final states in the reaction allows the 
calibration of the contribution of one- and two-body current contributions and makes it possible to assess which 
transitions are most sensitive to the presence of SRC with some confidence. This feature makes the "'^^O nucleus a 
prime target for this analysis, unlike the ^He nucleus which does not yield any bound states upon the removal of two 
protons. In Sec. II of this paper the essential ingredients of the description of the (e,e'pp) reaction and the calculation 
of the two-particle spectral function are summarized. The results are discussed in Sec. Ill, while conclusions are 
drawn in Sec. IV. 

II. CALCULATION OF THE (e,e'pp) CROSS SECTION 
A. Reaction mechanism 

The triple coincidence cross section for the reaction induced by an electron, with momentum Pq, where two nucleons, 
with momenta p'l and P2, are ejected from a nucleus is given, in the one-photon exchange approximation, by the 
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contraction between a lepton and a hadron tensor. If the effect of the nuclear Coulomb field on the incident and the 
outgoing electrons is neglected, the Lorentz condition for the MoUer potential and the continuity equation for the 
hadronic current make it possible to separate the longitudinal and transverse components of the interaction and to 
write the cross section as a linear combination of independent structure functions. For an unpolarized electron, after 
integration over the energy of one of the emitted nucleons {E2), the cross section is expressed in terms of six structure 
functions as |p|-pl]], 

d^cr 7re^ 

— Fv ilf /roc [ 2eL/oo + /ii - e (/i-i cos 2a + sin 2a) 



dp'Qd%dE[dn[dn'2 2q 

+ (1 + e) (/oi cos a + /oi sin a)] , (1) 

where e^/47r ~ 1/137, Pq is the momentum of the scattered electron, a is the angle between the plane of the electrons 
and the plane containing the momentum transfer q and p'l . The quantity 

1-^tan^-j (2) 

measures the polarization of the virtual photon exchanged by the electron scattered at an angle 9 and 

eL = -%e, (3) 
Q 

where — uP' ~ , with to — — and q — — p'^-, is the four- momentum transfer. The factor 

is the flux of virtual photons, fif — p'lE^p^Ei^ is the phase-space factor and 

Eb p'2 

is the inverse of the recoil factor. The quantity Eb is the total relativistic energy of the residual nucleus with 
momentum Pb = Q P'l P'2- 

The structure functions fw represent the response of the nucleus to the longitudinal (A = 0) and transverse 
(A — ±1) components of the electromagnetic interaction and only depend on lu, q, p[ p'2 and the angles 71 between q 
and p'l , 72 between q and p'^ and 712 between p'^ and p'2 ^ ■ They result from suitable combinations of the components 
of the hadron tensor [pljlll] and are thus given by bilinear combinations of the Fourier transforms of the transition 
matrix elements of the nuclear charge-current density operator taken between initial and final nuclear states 

J'^(q) = /(v|/f|>(r)|vI/i)e"'Mr. (6) 



These integrals represent the basic ingredients of the calculation. 

If the residual nucleus is left in a discrete eigenstate of its Hamiltonian, i.e. for an exclusive process, and under the 
assumption of a direct knock-out mechanism, the matrix elements of Eq. can be written as [p|jl2| 

J^{q) = J ipf{ri(Ti,r2(T2)J''{r,rieTi,r2cr2)i^i{rieTi,r2(T2)e"^"^drdridr2deTideT2. (7) 

Eq. contains three main ingredients: the final-state wave function ipf, the nuclear current and the two-nucleon 
overlap integral ipi . The bound and scattering states ipi and ip{ are consistently derived from an energy-dependent non- 
hermitian Feshbach-type Hamiltonian for the considered final state of the residual nucleus. They are eigenfunctions 
of this Hamiltonian at negative and positive energy eigenvalues, respectively [^,|ll[ . 

The same theoretical model for the exclusive (e,e'pp) reaction as in Ref. [12| is used, but here an improved treatment 
of the nuclear current and of the two-nucleon overlap integral has been adopted, as described below. 

In the final-state wave function ip{ each of the outgoing nucleons interacts with the residual nucleus while the 
mutual interaction between the two outgoing nucleons is neglected. The scattering state is thus written as the 
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product of two uncoupled single-particle distorted wave functions, eigenfunctions of a complex phenomenological 
optical potential which contains a central, a Coulomb and a spin-orbit term. The effects of an isospin-dependent 
term, to account for charge-exchange final-state interactions, were evaluated for the (e,e'pp) reaction in Ref. but 
negligible contributions were obtained in all the situations of practical interest. Thus this term is neglected here. 

The nuclear current is the sum of a one-body and a two-body part. The one-body part contains a Coulomb, a 
convective and a spin term. The two-body component is derived from the effective Lagrangian of Rcf. , performing 
a nonrelativistic reduction of the lowest order Feynman diagrams with one-pion exchange. In this approximation only 
processes with A-isobar configurations in the intermediate state contribute to the (e,e'pp) reaction. They produce a 
completely transverse current, J^. The operator form of was derived in Ref. [Q. It results from the sum of the 
contributions due to two types of processes, corresponding to the excitation and de-excitation part of the current. 
In the former case, the A is excited by photon absorption and then de-excited by pion exchange. The latter process 
describes the time interchange of the two steps, i.e., first excitation of a virtual A by pion exchange in a NN collision 
and subsequent de-excitation by photon absorption. For a pp pair they give |14 



J%jj{q, <Ti, <T2) = i 72rf ^ (2ifc T fc x a^^A x q Ga(V^) ^i'iD^ 



(8) 



where k is the momentum of the exchanged pion, m is the pion mass and the factor 7 collects various coupling 
constants, 7 = /■yNAfTrNNf-KNA/'m^- The dipole form factor 
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takes into account the electromagnetic form factor of the isobar, which corresponds to the isovector form factor 



M 



used in the static quark model |E5[. The propagator of the resonance, Ga, depends on the invariant energy ^/s of the 
A, which is different for parts I and II. For the de-excitation current ^/sJJ is approximated by the nucleon mass M 
and 



Ga = (Ma-M)-i, 
where Ma = 1232 MeV. For the excitation current we use |^6| 

M, 



where ^^snn is the experimentally measured invariant energy of the two outgoing nucleons. This gives 



(10) 



(11) 



Ga{V^) = Ma - - -rA(V^) 



(12) 



where the decay width of the A, Fa, has been taken in the calculations according to the parameterization of Ref. p7[ . 

The two-nucleon overlap integral tpi in (|^) contains the information on nuclear structure. For a discrete final state 
of the ^^C nucleus, with angular momentum quantum numbers JAf, the relevant part may be expressed in terms of 
relative and center-of-mass (CM) wave functions as 



V'i(ri(Ti,r2(T2) 



E 

nlSjNL 



CnlSjNL (l>nlSj{r)RNL{R) ^fsi^r, Tl, (T2) FL(^^fl) 



JM 



where 



ri - r2 



V2 



R 



V2 



(13) 



(14) 



are the relative and CM variables. Note that we follow the convention to denote lower case to relative and upper 
case to CM coordinate quantum numbers. The brackets in ( p^ indicate angular momentum coupling of the angular 
and spin wave function S of relative motion with the spherical harmonic of the CM coordinate to the total angular 
momentum quantum numbers JM. The CM radial wave function R is that of a harmonic oscillator [ p8[ , but the 
radial wave function 6 of relative motion includes a defect funtion in order to account for SRC flSll 



(kniSjir) = Rni{r) + Disj{r). 



(15) 
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These defect wave functions were obtained by solving the Bethe-Goldstone equation in momentum space for ^^O p9[ . 
For the present application these defect functions were Fourier Bessel transformed into coordinate space. This is not 
an exact procedure; the solution of the Bethe-Goldstone equation yields a non-local correlation operator which cannot 
strictly be represented by a local correlation function Disj of the form displayed in Eq. ([l5|). However, for the ^Sq 
wave, which is decoupled from other partial waves, the approximation is quite satisfactory. For the higher partial 
waves of the pp wave function the effect of SRC is relatively small due to the presence of centrifugal terms. 

The defect functions for the ^5*0 partial wave are displayed in Fig. |] for the Bonn-A, Bonn-C and Reid Soft Core 
potentials. One of the objectives of the present study is to investigate to what extent the differences between these 
defect functions are reflected in the calculated cross sections. 

The coefficients c in Eq. jl^ ) contain contributions from a shell-model space which includes the Is up to the 2plf 
shells. The framework within which this is done is basically the same as the one adopted in a recent calculation of 
the two-proton removal spectral function in momentum space JlSf . The main ingredients of this method are briefly 
presented in the next subsection. 
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FIG. 1. Defect functions (see Eq. (15)), multiplied by p = \p-^ — P2\/2, calculated for the ^5*0 partial wave by solving the 
Bethe-Goldstone equation in ^®0, by the method of Ref. [28]. Results are plotted for the Bonn-A, Bonn-C and Reid Soft Core 
potential. 



B. Structure amplitudes 



The guiding principle followed in the calculation of the structure amplitudes, which was presented earlier in Ref. 
I p^ , is the attempt of treating long-range and short-range correlations in a separate but consistent way. The effects 
of long-range correlations are determined by performing a nuclear structure calculation within a shell-model space 
including single-particle states which range from the Is up to the l/2p shell. Limiting the expansion in Eq. (|l3|) to 
this model space should be sufficient to accomodate long-range correlations. The effects of the strong short-range 
components of a realistic NN interaction, which would scatter the interacting nucleons into much higher shells, are 
taken into account by solving the Bethe-Goldstone equation using a Pauli operator which considers only configurations 
outside this model space. The solution of this Bethe-Goldstone equation yields the residual interaction of the nucleons 
inside the model space as well as the defect functions employed in Eq. (^5|) . The depletion of filled orbits by SRC is 
also incorporated in the shell-model space calculation by the energy dependence of the G-matrix interaction, which 
yields an energy dependent Hartree-Fock term in the self-energy [l9| ]. The fragmentation of one-nucleon removal 
strength is described by two-particle-one- hole and two-hole-one-particle terms in the self-energy S* in Tamm-Dancoff 
approximation p9| , |30[ |, with which the Dyson equation for the one-body propagator 

5a/3(w) = glpiuj) + gl^iuj)^s{uj)gspiuj) (16) 

7(5 

is solved. In Ref. [|l9| these dressed propagators were used to calculate the one-nucleon removal spectroscopic factors 
for the low-energy final states in -'^^N. The comparison with the results of one-nucleon knock-out experiments is then 
a first test of the quality of this ingredient in the calculation of two-nucleon removal amplitudes. The latter are 
contained in the Lehmann representation of the two-nucleon propagator G^^ 
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GiUA^) = (17) 
^(^^||(a^aa.)j||^]/}--^+^)(vI/'j-^+^||(atat),||vl/^) 



E 



uj — [h J — h ' ) + IT] 



UJ — (L — l^J ) — T-V 



\ ^ ^abJ^cdJ \ " ^cdJ^abJ 



/ ^ / r,n,A+2 TjiO.Ax , • / ^ z zpO,^ 7-nm,A-2\ 



E 



The symbols (. . . 1 1 . . . 1 1 . . .) represent the reduced matrix elements [pl|^3[ of the two-nucleon removal and addition 
tensor operators that are constructed by the angular momentum coupling of two one-nucleon addition and removal 
tensors aj^ and a^, where — (— )^°~™"a^Q, is the time reverse of a; —a denotes ^q, Jq, — ma} and a denotes 
basis states without the magnetic quantum number: a = {ua, la,ja}- 

The two-nucleon propagator is obtained by solving, within the shell-model space, the Bethe-Salpeter equation 
p3| , ^ for the two-nucleon propagator G^^ 

^1^375(^17*2,^3,^4)= (18) 

[9ati,{h -t'l)9t3y{h - ^2)] r^^KA(*i'*2'*3'*4)G'KA75(*3' *3, 4) , 

where T denotes the irreducible effective particle-particle interaction, which is here approximated by the G-matrix 
interaction which contains only propagation of particles outside the chosen model space. 

In the calculation of Ref. the spectroscopic factor for the removal of one nucleon from the p shell of ^^O turned 
out to be reduced by a factor 0.75 as compared with the independent-particle shell model. This is still about 10% 
larger than the factor 0.65 deduced from experiments [Q. We decided not to replace the calculated spectroscopic 
factor by the experimental ones in the dressed propagators. This means that the iwo-nucleon removal amplitudes 
that we obtain in the RPA with these dressed propagators |Q may be too large as well. This observation applies 
mostly to the noninteracting part of the two-particle spectral function represented by the first contribution to the 
two-nucleon propagator in Eq. ( [T^ ) . This term also yields a spurious contribution to the cross section for the one-body 
current contributions at small momenta jl^. The issue of interest here involves the effect of SRC which appear at 
higher momenta and the problem of spuriosity is not important. The overestimate may be much less severe for the 
interacting part of the spectral function (second term in Eq. (|l8|)) which yields the genuine SRC contribution to the 
cross section. In addition, such a factor, representing this overestimate, will be roughly the same for all the low-energy 
amplitudes involving removal of two protons from the p shell and therefore this uncertainty cancels in the comparison 
of relative magnitudes of amplitudes and cross sections for the low-energy final states in ^"^C. 

The shell-model two-proton removal amplitudes are expanded in terms of relative and CM wave functions for 
the initial state of the knocked-out pair. Summation over the contribution of the various configurations yields the 
coefficients c in Eq. (|l3|): 

^nlSjNL — 



Y.Y.^-)''^^+'+%2\ + l)'jSjah\sa sl s\{nlNL\\nJanbh\)\g j H , (19) 



with the notation j — y'2j+T and the nine-j and six-j symbols come from the angular momentum recouplings involved 

. . p . . . 

The most important amplitudes are listed in Table |. It is instructive to note that for these low-lying positive parity 

states the relative ^Sq wave is combined with a CM L = (for 0"*") or L = 2 (for 2+) wave, while the relative waves 
occur always combined with a L = 1 CM wave function. This was the basis of the global analysis of the experimental 
cross section in terms of ^5*0 and '^P removal contributions in Ref. The amplitudes for the 0+ states are presented 
at some length to illustrate the importance of the pairing interaction which mixes the shell-model configurations. 
Without this interaction, the lowest state would just correspond to the removal of two (dressed) protons from the pi 



6 



shell and the excited 0+ state to the removal from the ps shell. In that case the ^Sq removal cross section would be 
twice as large for the excited state as for the ground state. Due to the residual interaction the ground state becomes 
the strongest for removal, not only due to the coherent superposition of the p shell configurations but also the 
deep Is shell and the higher sd and pf major shells contribute. The contribution from these higher shells is much 
smaller for the 2+ states and completely negligible for the 1+ state. 

Another point to be mentioned is that the squares of the amplitudes add up to only about 0.6, as to be expected 
on the basis of the products of two one-nucleon removal spectroscopic factors (0.75)^. 
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FIG. 2. The differential cross section of the ^^0(c,c'pp) reaction as a function of the angle 72 for the transitions to the 
low-lying states in "C: 0^ (Sam = 22.33 MeV), 0+ (£2m = 32.08 MeV), 2+ {Eim = 30 MeV), 2+ (Bsm = 35.47 MeV), 1+ 
(£'2m = 33.64 McV). Eo = 584 McV, cu = 212 MoV, q = 300 McV/c, T{ = 137 McV and 71 = -30°. The defect functions for 
the Bonn-A potential and the optical potential of Ref. [37] arc used. Separate contributions of different relative partial waves 
are drawn. The contribution of the partial wave is very small for the 0"^ states and omitted from the figure. The solid 
lines give the total cross sections resulting from the contributions of all the relative states. 

III. TWO-PROTON KNOCK-OUT CROSS SECTIONS 

A. Relative magnitude of the contributions from one-body and two-body currents 

Of major interest in the (e,e'pp) studies is the question whether one may clearly identify the contributions from one- 
body and two-body currents and thereby study them separately. The part involving the one-body current is expected 
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to provide then an opportunity to probe SRC. These SRC, induced by the repulsive NN interaction, with a range 
of typicaUy 0.5 fm, wiU strongly affect the relative ^Sq wave function, but the short-range repulsion will have only a 
minor impact on the higher partial waves. For this reason a first inspection of the experimental data from NIKHEF 
has been made in Rcf. to estimate the relative contribution of ^5*0 and '^P pair knock-out in the cross sections for 
the lowest states of ^"^C. This estimate was based on the comparison of a simple factorization approximation of the 
cross section with the observed distribution of CM momenta of the knocked-out pairs (see also Ref. [0). Here we 
present the separate contributions of the ^5*0, ^Pj and ^Z?2 relative partial waves to the "'^^0(e,e'pp) cross sections 
for the low- lying states in ^*C. They are displayed in Fig. || for a specific kinematical setting that is included in the 
aforementioned NIKHEF data, with Eq = 584 MeV, 9 = 26.5°, uj = 212 MeV, and q ^ 300 MeV/c. The kinetic 
energy of the first outgoing proton T{ is 137 MeV. The missing energy i?2m = ^ ~ T[ — T2 — T^, where and Tg 
are the kinetic energies of the second outgoing proton and of the residual nucleus, respectively, has been taken in the 
calculations, for each transition, from a comparison with the experimental spectrum of ^^C psf but for the 2^ state, 
unidentified in the experimental spectrum, from the calculation of Ref. ||T^. The angle 71 is —30°, on the opposite 
side of the outgoing electron with respect to q. Changing the angle 72 on the other side the cross section can be 
explored at different values of the recoil momentum pg . The relationship between 72 and is shown in Fig. || for the 
transition to the ground state of ^"^C. Only small differences are obtained for the other states, owing to the different 
value of the missing energy. 
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FIG. 3. The recoil momentum pB as a function of 72 in the same kinematics as in Fig. g. 

In a factorized approach, where final-state interaction is neglected, pg is opposite to the total momentum of the 
initial nucleon pair. Thus in this approach the shape of the recoil momentum distribution is determined by the 
CM orbital angular momentum L of the knocked-out pair. This feature is not spoiled by final state interaction, 
which modifies the pair momentum. In fact in Fig. |^ the shapes of the angular distributions for different transitions 
and separate contributions of different relative states are determined by the corresponding values of L, indicated 
in Table |. The shape of the total result is driven by the component which gives the major contribution. Due to 
final-state interaction there is interference of different partial waves in the total cross section. In some cases it can be 
important, but in certain kinematical regions this is of minor importance, because there either one is much stronger 
than the other. 

The figures show that the cross section for the 0+ ground state, for the Oj and to a lesser extent also for the 2'^ 
state of ^''C, receive a major contribution from the ^Sq knock-out, as opposed to the higher lying states 1+, where 
only ^P waves contribute, and 2^, where the '^P waves are more prominent. This feature is in agreement with the 
experimental findings of Ref. The defect functions used in the calculations of Fig. ^ were those of the Bonn- A 
potential |lj]. The results for the Reid Soft Core potential have a similar qualitatitive behaviour for this case and 
therefore are not presented here. Calculations with the Bonn-C potential have not been performed, but from the 
shape of the defect functions shown in Fig. |l| we do not expect any significant difference with respect to the results 
obtained with the Bonn-A potential. 
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FIG. 4. The difFerential cross section of the 



^0(e,e'pp) reaction as a function of 72 for the transitions to the , ^2 



states in ^'*C in the same kinematics as in Fig. |^. Defect functions and optical potential as in Fie. ^. 
the one-body and of the two-body A-current are shown. The solid lines are the same as in Fig. 



OJ and 1+ 
Separate contributions of 



As already mentioned, one may hope that the one-body current and thus correlations yield the dominant contribu- 
tion to the cross section in some kinematical regions when the knocked-out pair is in a ^S'o state. The knock-out of 
and higher partial waves will proceed mainly through the two-body A-current. To illustrate to what extent our 
calculations support these expectations, we have plotted in Figs. and ||the separate contributions from the one-body 
and two-body current to the same total cross section as in Fig. For the 0+ states the contribution of the one-body 
current is much larger than that of the two-body current and the angular dependence has the s-wave shape typical of 
the contribution for these states. The results with the Reid defect functions have a similar shape but are a factor 
two smaller . In fact the range of relative momenta pj.e( = |Pi ^ P2I/2 probed in this region is w 1.5 fm~^, where the 
ratio of the Bonn- A and Reid ^6*0 defect functions is « 1.4, which gives a factor two in the cross section. For larger 
values of the recoil momentum the s wave becomes smaller, while the p wave becomes relatively more important. In 
the range of angles between 100° and 140°, where the recoil momentum is small, one may therefore probe correlations 
in the relative ^Sq wave function. 
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FIG. 5. The difFerential cross section of the ^®0(e,e'pp) reaction as a function of 72 for the transitions to the 2^ and 2j 
states in in the same kinematics as in Fig. |2[ Separate contributions of the one-body and of the two-body A-current are 
shown for the defect functions calculated with the Bonn-A and Reid potentials. The solid lines give the total cross sections 
resulting from the sum of the one-body and of the two-body A-current. Line convention and optical potential as in Fig. W. 



In sharp contrast to the 0+ states is the situation for the 1+ state. It is only reached by the knock-out of 
pairs and, as expected, the two-body current gives here by far the dominant contribution to the cross section. It will 
therefore be interesting to identify this cross section for the l"*", which is known to be at 10.3 MeV excitation energy. 

For the 2^ state we find that the one-body current gives a larger contribution than the two-body current, as 
opposed to the situation for the 2t state. This may be traced back to the large contribution of the ^5*0 partial wave 
for the 2^, as was shown in Fig. 0. For the 2 J especially dominates. However, the predicted dominance of the 
one-body contribution to the 2^ cross section depends on the defect functions used. This is shown by the comparison 
between the results obtained with defect functions from the Bonn-A and from the Reid potential in Fig. |5[ With the 
Reid defect functions the one-body current contribution is almost a factor two smaller than for the Bonn-A defect 
functions. This is not a general statement, but it turns out to be the case for the present kinematics. The cross section 
calculated with the two-body current is, as expected, only slightly affected by the choice of the defect functions. With 
the Reid defect functions the amplitudes from one- and two-body currents become of about the same size for the 2f 
state and the shape of the total cross section is determined by the interference of the two contributions. A similar 
result is obtained with the Bonn-A defect functions for the 2^ state. 
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FIG. 6. The differential cross section of the ^^0{e,e'pp) reaction as a function of 72 for the transitions to the 0^ and 2^ 
states in '^^C in the same kinematics and with same line convention as in Fig. ^ Defect functions and optical potential as in 
Fig. ^. The thin, mostly lower lines are calculated with the one-body current only. 

Next, we show explicitly how the amplitudes for knock-out from ^5*0 and higher partial waves are influenced by the 
A-current. This is plotted in Fig. ^ for the 0^ and 2^ states. The figures illustrate that indeed the ^5*0 knock-out 
amplitude is relatively little affected by the inclusion of the A-current, while this two-body current is a major factor in 
the knock-out of '^P and waves. This is a general result that has been obtained also in other kinematical situations. 
It can be understood if we consider the different role of the excitation and de-excitation part of the A-current. The 
excitation current, which has the energy-dependent A-propagator of Eq. (jl^), gives for energy transfer above 150 
MeV the dominant contribution of the A-current on "^P and waves. On a ^Sq pp pair the contribution of the 
excitation current is strongly reduced [^6| and becomes in our calculation of about the same size or even smaller than 
that of the de-excitation current, which is generally small. Thus the contribution of the A-current is generally small 
on a ^Sq pp pair, while it is dominant on '^P and pp pairs. The contribution of the waves to the total cross 
section is generally very small. So the proper place to study the two-body A-current in the (e,e'pp) reaction is where 
the ^P knock-out dominates, as in the 1+ and 2^ states, while SRC should be studied in the lowest states, where ^Sq 
knock-out dominates. Whether indeed one of these is dominant can be verified by inspection of the pair momentum 
distribution, as was illustrated in Ref. 111. 
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FIG. 7. The difFerential cross section of the ^^0(e,e'pp) reaction as a function of 72 for the transitions to the 0^ and 2^ 
states in "C, now in the same kinematics as in Ref. [12]: Eq = 475 MeV, u = 212 MeV, q = 268 MeV/c, Ti = 68 MeV and 
71 = 79.2°. Line convention as in Fig. kl 



B. Dependence on the NN potential and on the probed range of momenta 

In the discussion of Fig. || it was already indicated that especially the cross sections due to correlations and the 
one-body current are sensitive to the defect functions and thereby to the NN potential from which these were derived. 
For the range of relative momenta probed in the cross sections of Figs. |[^, the ^5*0 defect function of the Bonn- A 
potential is larger than that of Reid. In different kinematical situations it may be just the opposite. This appears to 
be the case for instance in the kinematics of Ref. ||l^ . In Fig. ^ it is shown that with that kinematics the contribution 
of the one-body current to the cross section for 0^ and 2^ is for most angles larger for Reid than for Bonn-A. The 
range of relative momenta probed here is on the average higher than in Figs. Around 7 w 65° the O]*" cross section 
is probed with prei « 2.1 fm~^. For the 2'^^ state the maximum around 7 ~ 90° corresponds to Prei ~ 2.2 fm~^. 

For really high relative momenta, above Prei ~ 3 fm^^, the contribution of the one-body current to the cross 
section will become systematically about a factor two larger for Reid than for the Bonn potentials. This is clear from 
the momentum dependence of the ^Sq defect wave functions that were shown in Fig. |l|. These might be probed in 
future experiments at TJNAF. Another possibility to discriminate between these potentials could be provided by the 
separation of structure functions. We discuss an example of this in the next subsection. 
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FIG. 8. The differential cross section of tlie ^®0(e,e'pp) reaction as a function of tlie recoil momentum pa for the transitions 
to the 0+, 2^ and 1+ states in "C in a super-parallel kinematics (71 = 0°, 72 = 180°) with Eo = 855 MeV, uj = 215 MeV and 
q = 315.89 MeV /c. The recoil-momentum distribution is obtained changing the kinetic energies of the outgoing protons. Line 
convention, optical potential and defect functions as in Fig. ^ Positive (negative) values of pa refer to situations where pg is 
parallel (antiparallel) to q. 



C. Separation of the structure functions /oo and /n in super-parallel kinematics 

The experimental separation of structure functions appears in general extremely complicated. The so-called super- 
parallel kinematics, where the knocked-out protons are detected parallel and anti-parallel to the transferred momentum 
q, is favored by the fact that only two structure functions, /oo and /n, contribute to the cross section, as in the inclusive 
electron scattering, and, as in that case, they can in principle be separated by a Rosenbluth plot Q. This kinematical 
setting has been realized in a recent experiment at Mainz | p5[ |. In this experiment, with an energy resolution of 
less than 1 MeV, different final states can be separated in the excitation-energy spectrum of the residual nucleus, in 
particular the 2+ states at 7.01 and 8.32 MeV. To compare the experimental results with our calculations, however, 
these two states should be considered as one state, the 2^, which is split up by the coupling to excitations of the ^^O 
core, that are very complicated and not included in our description. 

In Fig. ^ we display the cross sections for the Of ground state and the 2f and 1+ states in the super-parallel 
kinematics of the Mainz experiment, where 71 = 0°, 72 = 180°, Eq = 855 MeV, 9 = 18°, w = 215 MeV and q = 315.89 
MeV/c. The kinetic energy of the outgoing protons is changed in the calculations in order to explore different values of 
Pb- The figures show the decomposition into the different partial waves of the knocked-out pair. The recoil momentum 
distributions are similar to those shown in Fig. ^. The shapes of the different relative waves are determined by the 
corresponding value of L. The 0^ state is dominated for low values of pe, up to about 150 MeV/c, by ^So knock-out. 
At higher recoil momenta the contributions of ^5*0 and ^Pi knock-out become of the same order. We observe in this 
region that the total cross section may be lower than that given by the two separate contributions of ^5*0 and '^Pi 
states, owing to the negative interference of the two contributions. The 2^ state is dominated over the whole range 
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of recoil momenta by ^Sq knock-out, whose contribution is about a factor four larger than that of the other relative 
states. So in this kinematics the 2f seems to offer the best opportunity to study correlation effects. 

We do not display a decomposition into contributions from the one-body and two-body currents here, because the 
results are conceptually similar to those given in Figs. 4 and 5 and indicate the dominance of the one-body current 
for the 0+ and the 2^ states and of the A-current for the l"*" state. Moreover, the figures look quite similar to the 
ones shown here, i.e. the contribution of the one-body current is practically the same as that of the removal while 
higher partial waves come almost exclusively from the two-body current. This is illustrated explicitly for the 2f state 
in the last frame of Fig. ||. 

Essentially the same results as those shown in Fig. H, for the Bonn-A defect functions, are obtained with those of 
the Reid potential. In the latter case the one-body part is about 20% smaller, but otherwise the distributions are 
quite similar to those of Fig. ||. 
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FIG. 9. The structure functions foo and /n of the ^®0(e,e'pp) reaction as a function of pb for the transitions to the 
and 2^ states in ^"^C in the super-parallel kinematics of Fig. ^. Optical potential as in Fig. g. The solid and dashed lines are 
calculated with the defect functions of the Bonn-A and Reid potentials, respectively. 

A large difference between the results with the defect functions of Bonn-A and Reid potentials appears when a 
splitting into contributing structure functions /oo and /n is made. These results are plotted in Fig. ^. The transverse 
structure hmction /n appears to be insensitive to the defect functions. This is only partly due to the effect of the 
A-current. On the contrary the longitudinal structure function /oo, which is entirely due to the one-body current and 
thus to correlations, is much more sensitive than /n and the cross section to the choice of the defect functions. In 
this kinematics foo calculated with the Bonn-A defect functions is typically four times larger than calculated with the 
Reid defect functions. However, the experimental separation of the structure functions may be difficult, since /n is 
almost an order of magnitude larger than foo for the 0^ ground state. Also for the 1+, not shown in the figure, the /n 
structure function is found to be roughly five times foo with the Bonn-A defect functions and about twenty times foo 
with the Reid defect functions. Somewhat more favorable is the situation for the 2^ state, since here /n is only three 
times larger than foo at ps ~ 150 MeV/c, if the prediction with the Bonn-A defect functions turns out to be correct. 
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So this state may offer the best opportunity to determine the longitudinal structure function /oo experimentally. 

IV. SUMMARY AND CONCLUSIONS 

This work represents a combination of state of the art reaction description of the (e,e'pp) reaction together with 

a corresponding calculation of the two-micleon spectral function to produce results for cross sections mcasiired at 
NIKHEF and Mainz for the target. The description of the reaction includes both one- and two-body contribu- 
tions to the electromagnetic current. The treatment of final state interactions of the detected protons incorporates 
distortions (through an optical potential) for the individual particles but not their mutual interaction. Although the 
latter is expected to be unimportant for the cases of interest, this issue should be further studied in the future. The 
description of the two-body current involves a proper treatment of the dynamics of the intermediate excitation of the 
A-isobar before or after the absorption of the virtual photon. The two-nucleon spectral function (or two-nucleon over- 
lap function) has been obtained from a two-step procedure. The calculation of long-range correlations is performed in 
a shell-model space large enough to incorporate the corresponding collective features which influence the pair removal 
amplitude. The single-particle propagators used for this dressed Random Phase Approximation description of the 
two-particle propagator also include the effect of both long- and short-range correlations. In the second step that part 
of the pair removal amplitudes which describes the relative motion of the pair, is supplemented by defect functions 
obtained from the same G-matrix which is also used as the effective interaction in the RPA calculation. 

An important conclusion in this work concerns the predicted selectivity of the (e,e'pp) reaction involving discrete 
final states. Whereas the lowest 0+ and 2+ in ^"'C are predominantly reached by the removal of a pair other states 
at higher excitation energy mostly involve removal. The latter pair removal proceeds primarily via intermediate 
excitation of the A-isobar whereas the former is dominated by the one-body current mechanism. This feature is 
responsible for the calculated sensitivity in the cross sections to the treatment of short-range correlations where 
removal dominates. Short-range correlations induced by the Bonn- or Reid-potential may each yield larger cross 
sections than the other in certain kinematical domains. As a result, one will be able to study short-range correlations in 
this reaction successfully provided a sufficiently large set of kinematical conditions is explored including those available 
at TJNAF. The most promising extraction of the effect of short-range correlations shows up in the longitudinal 
striicture function which may be studied in the so-called super-parallel kinematics. Our study demonstrates that 
an intelligent choice of kinematics in exclusive (e,e'pp) experiments should allow the separation of the effects due 
to isobar currents and SRC for two nucleons with isospin T = 1. This success gives rise to the hope that a similar 
separation between two-body currents and SRC might also be possible in (e,e'pn) reactions. In this case one has to 
consider the competition between meson-exchange currents and SRC. The emission of a pn pair, however, probes the 
SRC for T = which are even stronger due to the presence of the tensor force. 
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TABLE I. Two-proton removal amplitudes from O for states of ^ C that are expected to be strongly populated in the 
^®0(e,e'pp) reaction. These are based on the Dressed RPA calculations described in Ref. [Q, within a model space of the 
Is up to the 2pl/ shells and with the G-matrix derived from the Bonn-C potential as an effective interaction. The quantum 
number p is the total number of harmonic oscillator quanta of the pair: p = 2n + Z + 2A'' + L (lower case for relative and upper 
case for center of mass motion) . For instance p = 4 indicates contributions from the sd shell and p = 6 from the pf shell. The 
energies of the listed states are largely known from experiments: 0'^ represents the ground state of ^*C, 2^ represents the sum 
of the 2+ states at 7.01 and 8.32 MeV Q, the 1+ is known ^ at 10.3 MeV. The 2 J was identified with a bump around 16 
MeV observed in Ref. The location of the Oj is less clear; the strength may be fragmented over several final states in the 
range between 12 and 14 MeV jl[. 
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